Isotropic-nematic phase equilibria in the Onsager theory of hard rods with length 

polydispersity 
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We analyse the effect of a continuous spread of particle lengths on the phase behavior of rodlike 
particles, using the Onsager theory of hard rods. Our aim is to establish whether "unusual" effects 
such as isotropic-nematic-nematic (I-N-N) phase separation can occur even for length distributions 
with a single peak. We focus on the onset of I-N coexistence. For a log-normal distribution we 
find that a finite upper cutoff on rod lengths is required to make this problem well-posed. The 
cloud curve, which tracks the density at the onset of I-N coexistence as a function of the width of 
the length distribution, exhibits a kink; this demonstrates that the phase diagram must contain a 
three-phase I-N-N region. Theoretical analysis shows that in the limit of large cutoff the cloud point 
density actually converges to zero, so that phase separation results at any nonzero density; this 
conclusion applies to all length distributions with fatter-than-exponentail tails. Finally we consider 
the case of a Schulz distribution, with its exponential tail. Surprisingly, even here the long rods 
(and hence the cutoff) can dominate the phase behaviour, and a kink in the cloud curve and I-N-N 
coexistence again result. Theory establishes that there is a nonzero threshold for the width of the 
length distribution above which these long rod effects occur, and shows that the cloud and shadow 
curves approach nonzero limits for large cutoff, both in good agreement with the numerical results. 



I. INTRODUCTION 

Rodlike particles such as Tobacco Mosaic Virus (TMV) 
in dilute suspension are known to exhibit a phase transi- 
tion with increasing density between an isotropic phase 
(I) with no orientational or translational order and a ne- 
matic phase (N) where the rods point on average along 
a preferred direction 1-5 . The main theoretical approach 
formulated to predict this phenomenon is the Onsager 
theory of hard rods 6 . Onsager assumed that the only in- 
teraction between the solute particles is of hard core type. 
The particles are modeled as perfectly rigid long, thin 
rods; non-rigidity as well as possible long-range attrac- 
tive potentials are neglected. Crucially, Onsager showed 
that the virial expansion truncated after the first non- 
trivial contribution becomes exact in the limit of long 
thin rods (the "Onsager limit"), i.e., for D/L — > where 
D is the diameter and L the length of the rods. The free 
energy then assumes a very simple form, because the sec- 
ond virial coefficient is just the excluded volume of two 
rods. The Onsager limit does however constrain the the- 
ory to low densities of order p ~ 0(1/DLq), and phases 
such as smcctics which occur at higher density cannot be 
predicted. 

In order to express the distribution of the non- 
conserved rod orientations, Onsager introduced the prob- 
ability P(fi) of finding a rod pointing along the direction 
O. Minimization of the free energy with respect to P(£l) 
results in a self-consistency equation for P(Q). Solving 
this in principle reduces the free energy to a function 
of the density p only, so that phase coexistences can be 
found by a standard double tangent construction. In or- 
der to avoid the complexity of the numerical solution 7 ' 8 
of the self-consistency equation, Onsager used a simple 
one-parameter variational trial form for P(Cl). Using this 
method Onsager 6 and, two years later Isihara 9 , were able 



to estimate the density at which the I-N phase transition 
occurs for different particle shapes. A similar approach 
was used by Odijk 10 , with a Gaussian trial function for 
P(f2). However, the numerically exact solution had by 
then already been obtained by Kayser & Raveche 11 . An 
alternative method, based on an expansion of the angular 
part of the excluded volume in terms of Legendre poly- 
nomials 11 , was used by Lekkerkerker et al. 12 . All of these 
approaches gave similar results for the properties of the 
coexisting isotropic and nematic phases. 

While being able to solve explicitly only the monodis- 
perse case, Onsager 6 already outlined the possible ex- 
tension of the theory to polydisperse systems, i.e., to 
mixtures of rods of different lengths and/or different di- 
ameters. Polydispersity has indeed been recognized as 
an important feature affecting experimental results 13,14 , 
and some attempts have been made to include it in the- 
oretical treatments 15-18 . A generic prediction is a pro- 
nounced broadening of the coexistence region with in- 
creasing polydispersity 19-22 , which is also observed ex- 
perimentally 14 . A second generic effect of polydispersity 
is fractionation, i.e., the presence of particles of differ- 
ent size in the coexisting phases; for rodlike particles, 
already Onsager 6 had predicted that the nematic phase 
would be enriched in the longer rods. Polydispersity can 
also result in more drastic and qualitative changes to the 
phase behavior, however. In particular, in systems with 
length polydispersity coexistence between two nematic 
phases (N-N) or one isotropic and two nematic phases 
(I-N-N) can occur. This has been observed experimen- 
tally 13 and predicted theoretically for bi- and tridispcrsc 
systems, i.e., mixtures of rods with two or three different 
lengths 16,17 . However, a detailed investigation of the ef- 
fects of full length polydispersity, i.e., of a continuous dis- 
tribution of rod lengths, on the Onsager theory remains 
an open problem. Perturbative approaches 20 ' 21 by their 
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nature cannot access qualitative changes to the phase 
diagram such as the occurrence of N-N or I-N-N coexis- 
tence. Our lead question for this paper is therefore: can 
N-N and I-N-N coexistence occur in length-polydisperse 
systems of thin hard rods? In cases where the length 
distribution has two or three strong peaks, one expects 
behavior similar to the bi- or tridisperse case, so that the 
answer should be yes. Much less clear is what to expect 
for unimodal length distributions, and this is the case 
that we will consider. 

We concentrate on the onset of isotropic-nematic phase 
coexistence coming from low density, i.e., on the isotropic 
cloud point; this can be calculated numerically with some 
effort using an algorithm which we have recently devel- 
oped 23 . Finding the phase behavior for higher densities 
inside the coexistence region would be substantially more 
difficult. We choose to start our analysis from a fat-tailed 
(log-normal) length distribution with a finite upper cut- 
off on rod length. This choice is inspired by the inter- 
esting results obtained by Sole 24 ' 25 for polydisperse ho- 
mopolymers, and by our recent investigation 26 of length 
polydispersity effects within the V2 Onsager model; the 
latter is obtained by a simplification of the angular de- 
pendence of the excluded volume of the Onsager theory. 
We showed that within this simplified model I-N-N coex- 
istence is indeed possible in a system with a log-normal 
(and hence unimodal) rod length distribution. The cloud 
curve, which gives the density where phase separation 
first occurs as a function of the width of the length dis- 
tribution, exhibits a kink where the system switches be- 
tween two different branches of I-N phase coexistence. 
The shadow curve, which similarly records the density of 
the incipient nematic "shadow" phase, has a correspond- 
ing discontinuity. Precisely at the kink in the cloud curve 
a single isotropic coexists with two different nematics, so 
that this kink forms the beginning of an I-N-N coexis- 
tence region. Both the cloud and the shadow curve were 
found to depend strongly on the rod length cutoff; in 
the limit of large cutoff the cloud and shadow curves ap- 
proach the same limiting form, which is universal for all 
length distributions with a fatter-than-exponential tail. 
The nematic shadow phase has rather peculiar proper- 
ties, being essentially identical to the coexisting isotropic 
except for an enrichment in the longest rods; the longer 
rods are also the only ones that have significant orienta- 
tional order. 

The above results for the V2 Onsager model suggest 
that also in the unapproximated Onsager theory a rod 
length distribution with a fat tail should have pronounced 
effects on the phase behavior. We will show numerically 
that the cloud curve indeed has a kink, and the shadow 
curve a corresponding discontinuity, demonstrating that 
the phase diagram contains a region of I-N-N coexistence. 
In fact, the effects of the fat-tailed length distribution 
are even stronger than for the Vi Onsager model, with 
the nematic shadow phase containing essentially only the 
very longest rods in the system. The numerical results 
leave open a number of questions, and we therefore sup- 



plement it with a theoretical analysis. We show that 
the assumption of a nematic shadow phase dominated 
by the longest rods is self-consistent, and are able to pre- 
dict that in the limit of large cutoff the density of the 
cloud point actually tends to zero: even though the av- 
erage rod length is finite, the presence of a tail of long 
rods drives the system to phase separate at any nonzero 
density. Motivated by these results, we finally revisit the 
case of rod length distributions with an exponential tail, 
using the Schulz distribution as an example. Numerical 
results show, surprisingly, that even here a regime occurs 
where the length cutoff matters and the nematic shadow 
phase contains predominantly the longest rods; a kink in 
the cloud curve again reveals the presence of an I-N-N 
coexistence region. (This is stark contrast to our results 
for the V2 Onsager model 27 , where the exponential tail of 
the distribution produces no unusual effects.) By return- 
ing to our theoretical analysis we find that the long-rod 
effects are weaker for the Schulz distribution than for the 
log-normal case: they only occur above a certain thresh- 
old value for the width of the rod length distribution, and 
the large-cutoff limits of the cloud and shadow densities 
above this threshold remain nonzero. 

The paper is structured as follows. In Sec. II we outline 
the extension to continuous length distributions of the 
Onsager theory and derive the phase coexistence equa- 
tions for the isotropic cloud point. Sec. Ill begins with 
a brief description of the numerical method we used to 
locate the cloud point; we then show our results for the 
phase behavior for log-normal length distributions with 
finite cutoff. In Sec. IV we outline our theory for the large 
cutoff limit and compare with some numerical results at 
finite large cutoff, finding good agreement. Finally, in 
Sec. V we turn to systems with a Schulz distribution of 
lengths, giving numerical results and sketching an appro- 
priately modified theoretical analysis which again makes 
predictions in good agreement with the numerics. Sec. VI 
contains a summary and a discussion of avenues for fu- 
ture work. In App. A we review in outline the high den- 
sity scaling theory for the monodisperse Onsager theory 
which we exploit in our analysis of the large cutoff limit . 
In App. B the main approximation underlying our the- 
ory for the log-normal distribution is justified, while the 
appropriate modifications for the Schulz distribution are 
sketched in App. C. 

II. THE POLYDISPERSE ONSAGER THEORY 

The Onsager theory with length polydispersity models 
a system of hard spherocylinders with equal diameters D 
but different lengths L. We introduce a reference length- 
scale Lq and write L = IL$ where I is a dimcnsionlcss 
normalized length. The Onsager limit is then taken by 
considering D /Lq —> at constant values for the normal- 
ized lengths I. From now on we will refer to I itself as 
the rod length unless stated otherwise; it can in principle 
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range over all values between and oo. 

The thermodynamic state of the system is described 
by the density distribution p(/,fi). This is defined such 
that p(l, ft) dl (dfl/Air) is the number density of rods with 
lengths in the range I . . .1 + dl and pointing along a di- 
rection within the solid angle dQ around SI. In terms of 
spherical coordinates, with the z-direction taken to be 
the nematic axis, we have d£l — sin 9 d6 dip and the den- 
sity distribution is independent of the azimuthal angle ip, 
p(l, fi) = p(l, 9). It can thus be decomposed according to 

p(i,9) = p(i)P(e\i) = P p(i)P(6\i) 

Here p(l) is the density distribution over lengths, 



- / dcos6 p(l, 



6) = / d0p(l,9) 



where we have introduced the shorthand 

d9 = —dcosO 
2 

The overall rod number density is 



p = J dlp{l) = J dld9p(l,0) 

so that P(l) = p(l)/p gives the normalized length distri- 
bution. From the definition of p(l) it also follows that the 
orientational distributions P(6\l) for rods of fixed length 
are normalized in the obvious way, J dO P(9\l) = 1. No- 
tice that the factor 4tt in the definition of p{l, f2) has been 
chosen so that for an isotropic phase one has the simple 
expressions P(9\l) = 1 and p(l) = p(l,9). 

We can now state the free energy density for the poly- 
dispersc Onsager theory (see Ref. 21). We use units such 
that k^T = 1 and make all densities dimcnsionless by 
multiplying with the unit volume Vq = (it/4:)DLq. The 
free energy density is then 

/ = / dlp(l) [In p(l)- 1] + / dldB p(l)P(9\l)\nP(9\l) 



dl dl' d9 d9' p(l)p(l')P(9\l)P(9'\l')ll'K(9, 9') (1) 



The first term gives the entropy of an ideal mixture, while 
the second term represents the orientational entropy of 
the rods. The third term is the appropriate average of the 
excluded volume (8/ 7r) Vb^'l SU17I (with Vo absorbed by 
our density scaling) of two rods at an angle 7 with each 
other. The kernel 6 ' 11 K(9, 9') results from the average of 
(8/tt)| sin7| over the azimuthal angles ip, ip' of the rods, 



K(9, 9' 



dip' dip 



sin 7 1 



— / — \/l — (cos 9 cos 9' + sin 9 sin 9' cos tpY 
7T Jo 2n v 



As in the monodisperse case, the orientational distribu- 
tions P{9\1) are obtained by minimization of the free en- 
ergy, Eq. (1); inserting Lagrange multipliers to enforce 
the normalization of the P(9\l), one finds 



P(9\l) 



Jd9> eW) 
ip(6) = - f dl'd9'p(l')P{9'\l')l'K{9,9') 



(2) 
(3) 



The conditions for phase equilibrium are that coexisting 
phases must have equal chemical potential fi(l) for all rod 
lengths Z, as well as equal osmotic pressure. The chemical 
potentials can be obtained by functional differentiation 
of the free energy (1) with respect to p(l). The orien- 
tational distributions P(9\l) do depend on p(l) but this 
dependence can be ignored because the P{9\1) are cho- 
sen to minimize /. Carrying out the differentiation and 
inserting Eq. (2) gives 



lip{9) - In / d9> e W fl '> 



Sf 
8 P (l) 

= hxp(l) + J d£P{9\l) 

dl'd9d9'p{l')P(9\l)P(9'\l')ll'K(9,9') (4) 
lnp(/) - In / d0e'*W (5) 



The osmotic pressure can be obtained from the Gibbs- 
Duhem relation, which for a polydisperse system reads 



11= [dl p(l)p(l) - f 



Inserting Eqs. (2) and (4) then yields 



U = p-- I dld9lp(l)P{9\l)ip{9) 



(0) 



A. Isotropic-nematic phase coexistence 

We now specialize the phase coexistence conditions 
to I-N coexistence, and then eventually to the isotropic 
cloud point, i.e., the onset of I-N coexistence coming from 
low densities. The isotropic phase will have 



^(010 = 1, <p I (6) = ip I 



-cip\, 



IT 



(7) 



where we have defined the first moment pi of the density 
distribution p(l, 9) 



Pi 



dllp{l) = p / dllP(l) = p(l) 



3 



which represents the scaled rod volume fraction. p\ = 
(Lq/D)4>. We have also used the fact that a uniform 
average of the kernel K(8,6') over one of its arguments 
is just an isotropic average over (8/tt)| sin , y|, giving 



dO' K(6, 6') = / d 1 sin 7 | sin 7 | 

n 2 Jo 

= -- = 2 = ci 

7T 4 



The equality of the chemical potentials (5) gives for the 
density distribution in the nematic phase 



p N (l)=p\l) dOe 1 ^ 



where 



(8) 



(9) 



The full density distribution over lengths and orienta- 
tions is therefore, using Eq. (2), 



p N (l,9) = p N (l)P N (9\l)=p I (l)e l9W 



(10) 



and the osmotic pressure (6) of the nematic phase can be 
rewritten as 



n- 



i 



dldO p\l)e la{e) - - / dld9lp l (l)e lg( - e) ^{9) 



(11) 



In the following we will concentrate on the isotropic 
cloud point, where the isotropic "cloud" phase starts 
to coexist with an infinitesimal amount of nematic 
"shadow" phase. At the cloud point the isotropic density 
distribution over lengths, p (I), therefore coincides with 
the overall density distribution of the system, p(°\l), 
which we call the parent distribution. The parent dis- 
tribution can be written as p l ' (0 = pP {0) {l), where 
p = fdlp(°'(l) is the overall parent number density 
and the normalized parent length distribution. 

Since all properties of the isotropic cloud phase are de- 
termined by the parent, we till drop the superscript "I" 
in the following. We will also take the parent distribu- 
tion to have average length (I) = 1; any other choice 
could be absorbed into the reference length Lq. This im- 
plies that the density and (scaled) volume fraction of the 
isotropic phase are equal, p\ = p. With this notation, 
the density distribution (10) of the nematic shadow is 
p N (l,6) = pPW(l)e l *W and fully determined by p and 
g(0). The function g(9) must obey 

9(0) = ~P I dlM' P i0) (l)e la ^lK(9,e') + c lP (12) 



as follows from Eq. (3) for ip N (0) together with Eq. (9). 
The cloud point density is the smallest value of p for 
which in addition the pressure equality is satisfied. Using 



Eq. (7) for the pressure of the isotropic and Eqs. (9) 
and (11) for that of the nematic, this condition reads 



p+\c lP 2 =p \ dld9P^(l)e 



f J dld6lpW(l)e ls W[g{6)-c lP } (13) 



B. Fat-tailed rod length distributions 

So far everything is general and applies to any par- 
ent length distribution P(°\l). Let us now focus on 
the case of a parent distribution with a fat, i.e., less 
than exponentially decaying tail for large I. At the 
cloud point we have from Eq. (8) the density distribu- 
tion p N (l) = pP (0) (l) Jd9e tgW in the nematic shadow 
phase, or if we isolate the value of g{9) at 9 = 



p N (l)=pP^(l)e l9W 



d6 



oW)-g(0)] 



(14) 



In a nematic, one expects g(9) < g(0) and therefore the 
angular integral reduces to a less than exponentially vary- 
ing function of I. Moreover, g(0) is expected to be pos- 
itive, since the nematic phase should contain the longer 
rods. The nematic density distribution p^(l) is therefore 
exponentially diverging for large I whenever the normal- 
ized parent distribution p(°) (I) decays less than exponen- 
tially. In order to ensure finite values for the density and 
volume fraction of the nematic phase, we thus need to 
impose a finite cutoff l m on the length distribution; un- 
less otherwise specified, all Z-intcgrals will therefore run 
over . . . lm from now on. The presence of a cutoff is 
of course also physically reasonable, since any real sys- 
tem contains a finite largest rod length. Nevertheless, we 
will later also consider the limit of infinite cutoff, which 
highlights the effects of the presence of long rods. 



III. NUMERICAL RESULTS FOR THE ONSET 
OF I-N COEXISTENCE 

A. Numerical method 

A numerical determination of the isotropic cloud point 
involves the solution of the two coupled equations (12,13) 
for p and g(9). In an outer loop we vary the density until 
the smallest p that satisfies the pressure equality (13) is 
found; we use a false position method 28 . The nontrivial 
part of the algorithm is the inner loop, i.e., the solution 
of the functional equation (12) for g(9) at given p. An 
iterative method inspired by the one used by Herzfeld 
et al. 7 for the monodisperse case turns out to converge 
too slowly in the presence of polydispersity. We therefore 
choose to represent g(9) by its values gi = g(9i) at a set 
of n discrete points 9^ the values of g(9) for 9 ^ 9i are 
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then assumed to be given by a cubic spline fit 29 through 
the points (6i,gi). This turns the functional Eq. (12) 
into a set of n nonlinear coupled equations which can be 
solved by e.g. a Newton-Raphson algorithm 28 . To keep 
n manageably small while keeping the spline representa- 
tion accurate, a judicious choice of the 9i is important. 
We exploit the symmetry g(9) = g(jr — 9) and choose 
a nonlinear (geometric) spacing of the Oi over the range 
. . . 7r/2, with more points around the origin where g{9) 
is least smooth. 

The numerical method outlined above was also used 
in our recent investigation of the exact phase diagram 
for the bidisperse Onsager theoery 23 . The only modi- 
fication is in the /-integrations, which in the bidisperse 
case became simple sums over the two rod lengths. In 
principle the additional integration could make the cal- 
culation substantially slower; however, one notices that 
in Eqs. (12) and (13) only the two /-integrals 



h{x) 



dl P^\l)e lx 



and 



h\x) = dl P (0 \l)lt 



are needed. We therefore precompute these functions and 
store them once and for all as cubic spline fits which can 
be evaluated very efficiently. 

As an alternative to the approach above we also consid- 
ered calculating g(9) by minimizing an appropriate func- 
tional. If p l (l) and II 1 are the chemical potentials and 
osmotic pressure of the isotropic parent phase, then one 
easily sees that a local minimum of the functional 

Z[p(l,6)} = f[p(l,e)]- f dldfffi 1 (l)p(l,9)+Il I (15) 



corresponds to a phase p(l, 9) with the same chemical 
potentials as the parent. Inserting for p(l,0) the known 
form (10) for the nematic density distribution, S turns 
into a functional of g{9), and the condition for a local 
minimum becomes equivalent to Eq. (12). S[g(#)] al- 
ways has g{9) = as a minimum, corresponding to the 
isotropic parent itself; but close to the onset of phase co- 
existence an additional nematic solution with g(9) ^ 
appears. Notice that, geometrically, S is a tilted version 
of the free energy for which the tangent (hyper-)plane at 
the parent is horizontal (3 = 0). At phase coexistence, 
i.e., for a parent with the cloud point density, the tangent 
plane also touches the nematic and so the isotropic and 
the nematic minimum of S are both at "height" 5 = 0. 

Numerically, one could minimize S by again represent- 
ing g{9) as a spline through a finite number of points 
.9; = g{0i) and then minimizing the resulting function 
of the gi. In general this turns out to be no easier 
than the solution of the similarly discretized version of 
Eq. (12). However, the minimization approach is useful 



when g{9) assumes a simple parametric form. We will 
see later that this is indeed the case for large cutoff l m , 
with g{9) being well approximated by the two-parameter 
form g(9) = a — b sin 9. Inserting this into S and minimiz- 
ing over a and b then gives an approximate solution for 
g(9); when used as a starting point for g(9), this makes 
it significantly easier to converge the numerical solution 
of the discretized Eq. (12) described above. 



B. Results for log- normal length distribution 

With the numerical method described in the previous 
section, it is possible to solve for the onset of isotropic- 
ncmatic phase coexistence for in principle arbitrary par- 
ent length distribution. We choose here a specific fat- 
tailed length distribution, the log-normal, which has al- 
ready given interesting results in polymers 24 and in our 
previous analysis of the Vi Onsager model 26 . The log- 
normal distribution has the form 



P<°>(Z) = 



1 1 



V2ttw 2 I 



exp 



Qnl-fj,) 



2 1 



2w 2 



(16) 



with finite length cutoff l m . The quantity w, which tunes 
the width of the distribution, is fixed by the normalized 
standard deviation a (in the following referred to as poly- 
dispcrsity) 

_2 (I") (I)' 



(I) 2 

to be w 2 = ln(l + a 2 ). The second parameter p is deter- 
mined so that the parent has average length (I) = 1, giv- 
ing p = —w/2. Notice that with these choices, the parent 
length distribution is normalized and has the desired mo- 
ments (I) = 1 and (I 2 ) = 1+a 2 only in the limit of infinite 
cutoff l m . The deviations for finite cutoffs are small even 
for relatively modest l m , however. For instance, at cutoff 
l m = 50 and a = 0.5, the integrals J ' m dl PP(°)(Z) for 
n = 0,1,2 differ from their respective values at infinite 
cutoff, 1, 1 and 1.25, by values of order 10 -17 , 10~ 15 and 
10~ 13 . Since we will not consider smaller cutoff values 
below, these small corrections can safely be neglected. 

In our previous analysis of the V2 Onsager model 26 , 
we observed that for log-normal length distributions the 
cloud curve has a kink, and the shadow curve a corre- 
sponding discontinuity; at the kink, the isotropic phase 
is in coexistence with two distinct nematics, so that the 
phase diagram must contain a region of I-N-N coexis- 
tence. In the V-2 Onsager case, the simplicity of the 
model actually allowed us to compute the complete phase 
diagram and locate the three-phase I-N-N region explic- 
itly. For the full Onsager theory treated here we can only 
find the cloud and shadow curves at present, not the full 
phase diagram; nevertheless, a kink in the cloud curve 
will again imply the presence of an I-N-N region in the full 
phase diagram. In Fig. 1 we show the cloud and shadow 
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completely dominated by the longest rods in the parent 
distribution. In the rescaled volume fraction (pi) repre- 
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FIG. 1. (a) Number density p of the cloud phase for length 
cutoff l m — 50 (solid) and l m = 100 (dashed), plotted against 
the polydispersity a on the y-axis. Notice the kinks in the 
two curves, which imply the presence of a three-phase I-N-N 
coexistence region in the full phase diagram. The kinks cor- 
respond, as they should, to discontinuities in the shadow 
curves (b). Both cloud and shadow curves are strongly cut- 
off-dependent, moving towards lower densities as l m increases. 



curves, i.e., the number density of the isotropic cloud (a) 
and nematic shadow (b) plotted against polydispersity, 
for a log-normal parent with two different cutoffs. The 
presence of the kink in the cloud curves, and the cor- 
responding discontinuity in the shadow curves, is clear 
evidence of the presence of an I-N-N coexistence region 
starting at the kink of the cloud curve. The positions 
of the kink and discontinuity, respectively, as well as the 
shapes of the cloud and shadow curves above them, show 
a strong dependence on the cutoff length; both curves 
move to significantly smaller densities as l m increases. 
For polydispersities a below the kink, on the other hand, 
the number of long rods is too small to have a significant 
effect on the phase separation and one has essentially 
cutoff-independent behavior that connects smoothly with 
the monodispcrse limit a = 0. These observations are in 
qualitative accord with our earlier results for the Vi On- 
sager model with the same length distribution 26 . 

Moving across the discontinuity from below, the 
shadow curve jumps from a "normal" nematic branch to 
an unusual nematic phase which, as we will see below, is 




FIG. 2. Scaled volume fraction p ± of the nematic shadow 
phase at l m = 50 (solid) and l m = 100 (dashed). Notice that 
the discontinuity of the two curves is now much wider than 
in the number density representation of the shadow curves in 
Fig. 1(b). 



sentation of the shadow curve shown in Fig. 2 the differ- 
ent characteristics of the two nematic phases are clear. 
While at low polydispersity pf is of order unity, it jumps 
by two orders of magnitude on crossing the discontinu- 
ity. This shows that the fractionation effect which one 
normally encounters in polydisperse systems, with the 
long rods found preferentially in the nematic phase, be- 
comes extreme here. A plot of the average rod length 
in the nematic shadow against polydispersity (Fig. 3) 
in fact shows that above the discontinuity, the nematic 
phase contains almost exclusively rods of length close to 
the cutoff length. Since the longer rods are generically 
expected to be more strongly ordered, one should then 
also find that the nematic shadow phase has very strong 
orientational order. This is indeed the case: the ori- 
entational order parameter S, defined as an average in 
the nematic phase over the second Legendre polynomial 
V 2 {cos9) = (3cos 2 6>- f)/2, 

S = [ dl d9 P N (Z)P N (9\l)V 2 (cos 9) 



is almost indistinguishable from f above the discontinuity 
in the nematic shadow curve, implying that the typical 
angles 6 which the rods make with the nematic axis are 
very small. This also implies that the rods we consider 
must be rather thin for Onsager's second virial approxi- 
mation to be valid: for monodisperse rods 30 , the criterion 
is D/L <C 6. We return to this point in Sec. IV C. 

The results shown above for the finite cutoff regime 
leave open a number of questions. For example, we ob- 
served that both the isotropic cloud and nematic shadow 
curves move to lower densities as the cutoff increases, 
but the modest Z m -values used are too small to deter- 
mine whether the curves will converge to a nonzero limit 
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FIG. 3. Average length (Z) N = Pi/p N in the nematic 
shadow phase against polydispersity, for cutoffs l m = 50 
(solid) and l m = 100 (dashed). Above the discontinuity, 
fractionation becomes extreme. In fact, at the discontinu- 
ity the nematic phase appears to be composed only of rods 
with lengths of order the cutoff length. At the top of the 
plot, for very large polydispersities where the parent length 
distribution becomes somewhat more uniform, some shorter 
rods are also included in the nematic phase and reduce the 
average length. 



as l m grows large or instead approach zero. One would 
also like to ascertain whether the average rod length in 
the nematic shadow really tends to l m for large cutoffs, 
as suggested by Fig. 3, and what happens to the rescaled 
nematic volume fraction pf in the limit (Fig. 2 suggests 
that it might become large). In the next section we there- 
fore turn to a theoretical analysis of the limit l m — > oo, 
which will clarify all these points. 



IV. THEORY FOR FAT- TAILED 
DISTRIBUTIONS WITH LARGE CUTOFF 

Above we saw that, at the onset of isotropic-nematic 
phase coexistence in systems with log-normal length dis- 
tributions, the nematic shadow phase appears to be 
dominated by the longest rods in the distribution, with 
lengths I « l m . In this section we will construct a con- 
sistent theory for the phase behavior in the limit of large 
cutoffs l rn based on this hypothesis. From this we will 
be able to extract the limiting dependence on l m of the 
cloud point density p, the density p N and (scaled) vol- 
ume fraction p^ of the nematic shadow phase, and the 
function g(6) determining the oricntational ordering of 
the nematic phase. At the end of this section we will 
then compare these predictions with numerical results 
obtained for finite but large cutoff. 



A. Dominance of long rods in the nematic 

Recall that, in principle, we need to solve Eqs. (12) 
and (13) for g(6) and p to determine the isotropic 
cloud point and the properties of the coexisting nematic 
shadow phase. It will be useful to isolate, in Eq. (12), the 
contribution e' 9 ' ^ which determines the divergence of the 
nematic density distribution (14). Define the function 



h(9) = l m \g(0) - 9(0)] 



(17) 



with the sign chosen so that h(9) should be non- negative 
for all 9. Eq. (12) evaluated at 6 = then gives 

g(0) = -p j dl pWQ)e lB <-°H J dff e-^ l ^ h ^K(0,6') 
+ c lP (18) 
and subtracting from Eq. (12) yields 

h(6)=pl m f dl P(°\l)e l9 (°h 

x [d9> e' (l/lm)h(e ' ) [K(6, 9') - K(0, 9')] (19) 



Together, Eqs. (18) and (19) for g(0) and h(9) are of 
course equivalent to Eq. (12) for g(9). 

To formalize the assumption that the nematic phase is 
dominated by the longest rods, consider now the nematic 
density distribution (14) which in our new notation reads 



,p(o)(/) e fo(o) / d0 e -d/U)h(0) 



(20) 



If the exponential factor exp[Z<7(0)] is large enough for the 
nematic indeed to be dominated by the longest rods, we 
can replace the weakly varying (at most as a power law 
in /) angular integral by its value at I = l m , giving for 
the nematic density p N = J dl p N (l) 

dl p(°)(Z) e ^(°) J d9e~ h ^ (21) 

We can now make the same approximation in Eq. (18): 
the weakly varying factor in the ^-integral is I times the 
angular integral, and replacing this by its value at / = l m 
yields 

5(0) = -pl m J dl P (0 \l)e l9i0) J d£ e- h ^K(0,6)+c lP 

(22) 

In Eq. (19) for h(9), the analogous approximation gives 

W- A 

x / dff e-' l < e ') [K(e, 9') - K(0, 9 1 )] (23) 
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Finally, consider the expression on the r.h.s of Eq. (13) 
for the osmotic pressure IT N in the nematic. The first 
term is the ideal contribution, i.e., the nematic density 
p N , while inserting the definition (17) in the second term 
gives 

jjN N 



-h(0) + g(0) - c lP 



(24) 



Replacing I by l m in the weakly varying terms then yields 
n N = p N -f / d/P (O) (0e' s(O) 

x fd9e- h ^[-h(e) + l m (g(0)~c lP )} (25) 



We show in App. B that a posteriori, the approximation 
of dominance of the long rods can be justified in all cases 
above, with the contributions to the /-integrals from rod 
lengths I <C l m becoming negligible for l m — > oo. 



B. The large-cutoff scaling solution 



We have now got four equations to be solved for 
P N ; 5(0)i h(9) and p, appropriately simplified using 
the assumption that the nematic phase is dominated 
by the longest rods in the parent distribution. These 
are Eqs. (21), (22) and (23) and the pressure equality 
n N = p+(ci/2)p 2 with n N from Eq. (25). Using Eq. (21), 
Eq. (23) can be written as 

h{9) = p eS - — (26) 



f d e> 



-h(e> 



Here we have defined 



Pcff 



(27) 



which is just the dimensionless density of the nematic 
phase, with the factor arising from the fact that, 
since the nematic is effectively monodisperse with I = l m , 
one should use l m Lo rather than Lo in the definition of 
the unit volume Vq (see before Eq. (1)). In the form 
above, Eq. (26) is identical to Eq. (A4) in App. A for a 
monodisperse system at dimensionless density p c g. An- 
ticipating that p c ft will be large, an assumption to be 
checked a posteriori, we then deduce immediately that 
h(9) will be given by the high-density scaling solution 
sketched in App. A, h{&) = h(t), with the scaling vari- 
able t = p ff sin0. The scaling function is determined by 
Eq. (A8) 



h(t) = {K{t,t') - K{0,t')) t , 



(28) 



where the average over t' is over the normalized distribu- 
tion 



P(t) = -te 
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-h(t) 



7 = dt te 



-h(t) 



(29) 



and K{t,t') is the small-angle scaling form of the ker- 
nel K(9,9') given in Eq. (A7) in App. A. The range 
of all averages and integrals over t and t 1 can be taken 
as ... 00 (rather than . . . p c g) for large p e s since the 
largc-i regime gives only a negligible contribution. 

With the form of h(9) determined up to the single pa- 
rameter p e ff = /° N ^mi j us t three unknowns p, p N and g(0) 
now remain, to be determined from Eq. (21), Eq. (22) 
and the osmotic pressure equality. We now simplify 
these relations further, making use of the fact that in 
all angular integrals involving the factor exp[— h(9)] = 
exp[— h(p Q g sin 6)] only small angles 9 ~ l/p c ff contribute 
significantly. Physically, this means that the rods in the 
nematic shadow phase, which is at high dimensionless 
density p c ft, have strong orientational order. For such 
small 9 we can set sin 9 ~ 9 and transform everywhere to 
the scaling variable t = Pe ssm9 f=a p c g9. Eq. (21) then 
becomes, using the definition (29) 



p n = p [ di p(°)m e w°) / 
J J pis 

= £^L [ dlP^(l)e 1 ^ 

PeS J 



-h(t) 



(30) 



Eq. (22) can be similarly transformed and, using 
^(0,0) = (8/tt) sin6> = (8/ir)t/p eS , reads 



5(0) 



-Pi; 
f cip 



dl P(°)(0e i ^°) / ^*e-fcW (- — 



(31) 



Comparing with Eq. (30), and using p N Z m /peff = f/^m, 
this can be written in simpler form as an average over 
the distribution (29), 



5(0) 



7T lm 



Cip 



(32) 



Finally, the expression (25) for the osmotic pressure in 
the nematic can also be simplified by using that 9 is 
small and transforming to the scaling variable t. Insert- 
ing Eq. (32), this gives 



n- 



f dtt e _~ h[t) 


h(t) + -(t) t 


1 Pis 





= p* + t[ di p(°)(zys(°)4- 

2 J Peff 



(h(t)) t + -(t) t 
IT 



In App. A we show (Eq. (A15)) that the scaling proper- 
ties of the high-density limit imply that the constant in 
the square brackets has the value 4; using also Eq. (30) 
we then get the simple result 



n N = 3 P N 



(33) 
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This is identical to Eq. (A12) derived in App. A for 
monodisperse nematics at high (dimensionlcss) densities, 
and therefore consistent with our assumption that the ne- 
matic shadow phase behaves as an effectively monodis- 
perse system. The isotropic parent phase has pressure 
If 1 = p + (ci/2)p 2 , so that the pressure equality takes the 
form 



1 / , ci 



(34) 



We can now proceed to determine the ^-dependence 
of p, p N and g(0). Multiplying Eq. (30) by (p N ) 2 and 
inserting Eq. (34) gives 

The ^-integral will again be dominated by the longest 
rods, so that we can set P(°)(Z) = P(°)(I ro ) to leading 
order (see Rcf. 26) to obtain 



27 V 2 V 



J m g(0) 
5(0T 



Inserting Eq. (32) to eliminate <?(0), we finally get a non- 
linear equation relating p and l m 



27 7 V 2 H 



ciP , 

7T („ 



e W*)(t) 



To obtain the asymptotic solution for large l m , we antic- 
ipate that p will vary no stronger than a power law with 
l m ; for all fat-tailed parent distributions except those 
with power-law tails, the dominant Z m -depcndence on the 
r.h.s. will then be through the factor 1/P(°'(l m ). Taking 
logs we have 



l m c lP =-\nP^(l m ) + 0(\nl m ) 



(35) 



Specializing to log-normal parent distributions, with 
lnP( 0) (Z) = -(ln 2 /)/(2u) 2 ) to leading order, we finally 
get 



In 2 L 



P = 



2ciW 2 l r , 



O 



In l n 



(36) 



showing that p indeed varies with l m as a power-law (with 
logarithmic corrections). Our theory thus predicts that 
the isotropic cloud point density converges to zero for 
large cutoffs; in the extreme limit of a log-normal par- 
ent distribution with an infinite cutoff, phase separation 
would occur at any nonzero density. From Eq. (34), the 
density of the nematic shadow phase likewise vanishes for 
large cutoff, with p N = p/3 to leading order. For g(0) we 
have from Eq. (32) 



5(0) 



In 2 l„ 

2wH r 



O 



In L 



7T l m 



(37) 



and the last term of 0(l/l m ) is subdominant compared 
to both the first term and the first order correction 
O ((In l m )/l m ). We can also now write down the whole 
function g(9), using g(6) = g(0) - h(0)/l m = g(0) - 
h(p N l m sin 9) /l m . From App. A we know that the scal- 
ing function h is given by h(t) = (8/w)t up to correction 
terms of 0(1); once multiplied by l/l m , these give only 
subleading corrections to g(9). We thus find to leading 
order 

g(9) = g(0) - ^-h(p N l m sin 9) ~ a - bsin9 (38) 

where, using Eq. (36) and the leading-order relation p N = 
p/3 



a ee ,g(0) 



In 2 In 
2w 2 l r 



O 



\nL 



In 2 L 



b = -l m p N =- - 

it Sir Zciw^ 



0(lni 



(39) 
(40) 



The dominance of the long rods in the nematic phase thus 
results in a very simple form for g(9), with a vanishing 
and b slowly diverging in the limit l m — > oo. 

Having obtained the desired predictions from our the- 
ory, we can now also verify that the assumption of a 
large dimensionlcss density p Q ff for the nematic phase 
is justified. From Eq. (36) and the fact that to lead- 
ing order p N = p/3 one has p cS = p N l m ~ l m In 2 l m 
for a log-normal parent distribution, and this indeed be- 
comes arbitrarily large as the cutoff l m increases. The 
self-consistency of the other assumption which we made, 
i.e., the dominance of the long rods, is verified in App. B. 



C. Validity of Onsager theory 

Before comparing our theoretical predictions with nu- 
merical results at finite cutoff, we briefly assess the limit 
of validity of Onsagcr's second virial approximation. For 
a monodisperse system, an analysis of the scaling of the 
second and third virial coefficients 30 shows that in the 
nematic phase typical rod angles 9 with the nematic 
axis have to be 3> D/L, with D and L the diameter 
and length of the rods, for the truncation after the sec- 
ond virial contribution to be justified. In our situation, 
the nematic phase is effectively monodisperse with (un- 
normalizcd) rod length Lol m , so the condition becomes 
9 ^> D/(Lol m ). We showed above that the typical angles 
scale, for large cutoff l m , as 9 ~ l/p e s = l/(p N ^m)' so 
that the second virial approximation breaks down when 



D 

Loin 



„N/2 



(41) 



Now p N scales as (\n 2 l m )/l m , so this becomes In 2 l m ~ 
Lq/D. This shows that fairly large values of the aspect 
ratio Lq/D of the "reference rods" are necessary in order 
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for the theory to be valid for large cutoffs. For instance, 
one would need L /D ~ 50 for a cutoff l m = 1000. The 
longest rods are then very thin indeed, with L^l m /D ~ 
50,000. 

A more physically intuitive interpretation of the above 
condition is that it corresponds to the requirement of 
having a rod volume fraction <p <C 1. For monodisperse 
rods, we have (see App. A) that at large dimcnsionlcss 
density p, ~ p' 1 ~ (tfDN/V)- 1 . The limit of validity 
of the Onsager theory is therefore given by D/L ~ 9 <~ 
V/{L 2 DN) or 1 - {LD 2 N)/V ~ <f). The same is true 
for our calculation above: the volume fraction of the ne- 
matic phase is </> N = (D/L a )p^ ~ (D/L a )l m p N . The 
condition (41) thus again becomes </> N ~ 1, and we need 
N <C 1 for the second virial theory to be valid. 



g(0) = cip, while for smaller cutoffs deviations from the 
asymptotic theory appear as expected. 

For the parameter b, our theory predicted b = 
(8/37r)Z m p in the limit of large l m . In Fig. 5 we plot 
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D. Comparison with numerical results 

We now compare the theoretical predictions obtained 
above for the limit l m — > oo with numerical calculations 
for finite but large cutoff. Our numerical results will be 
able to confirm only the leading terms of the scaling solu- 
tion, since sub-leading corrections (e.g. to the result (38) 
for g(9)) can arise from the regime of very small angles 
9 ~ 1/peff = l/(p N ^m) which we cannot resolve numeri- 
cally. 

1.2 I 1 1 1 1 1 ^ 
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FIG. 4. Parametric plot of g(0) against p, for a log-normal 
distribution with a — 0.5 and a range of cutoffs between 50 
and 3000. At small p (large l m ) the numerical results (solid) 
are in very good agreement with the theoretically predicted 
asymptotic relation g(0) = cip (dashed). 

We begin by checking the predicted relations between 
the cloud point density p and the other quantities we 
have analysed theoretically, namely the parameters a and 
b specifying the leading behavior (38) of g(9), and the 
nematic shadow density p N . In Fig. 4 we plot <?(0) = a 
against p for a range of cutoffs l m between 50 and 3000. 
At large cutoff, i.e., at small p, the numerical results 
are clearly seen to approach the theoretical prediction 



20 ^ 1 1 1 1 1 1 

20 30 40 50 60 70 80 

(8/37r)/ m p 

FIG. 5. Parametric plot of b against the theoretically pre- 
dicted value (8/37r)Z m p for a log-normal distribution with 
o — 0.5 and cutoff between 50 and 3000. The convergence 
of the numerical results (solid) to the theoretical prediction 
(dashed) for large cutoff, i.e., large l m p, is clear. 

the numerically obtained b against this theoretical pre- 
diction, for a range of cutoffs between 50 and 3000. At 
large l m p the convergence to the theoretical solution is 
clear, while at finite cutoff (small l m p) the value of b lies 
rather below the theoretical expectation. This is not sur- 
prising, as b is expected to diverge with l m like In 2 l m and 
therefore rather slowly; at finite cutoff, then, correction 
terms will be rather important. 

Based on the fact that the nematic pressure obeys the 
simple relation II N = 3p N in the large cutoff limit, our 
theory also predicts that p N should be related to p by 
p N ='[p+(ci/2) j o 2 ]/3= (p + p 2 )/3. A plot ofp N against 
p for the same range of cutoffs as above (Fig. 6) clearly 
shows that this relation is satisfied in the limit of large 
l m , i.e., small p. In fact, deviations from the predicted 
scaling are rather small already at modest l m (here in the 
range 50 < l m < 3000). This shows that the dominance 
of the long rods, demonstrated also by the fact that the 
average length of the nematic phase is very close to l m 
(Fig. 3), appears quite early on. Even at l m = 50, not 
only is the nematic phase composed almost entirely of 
rods of order l m , but the rods are also sufficiently strongly 
ordered to make our scaling solution a good approxima- 
tion. 

Having confirmed the relations between a, b, p and p N , 
the last prediction to verify is the variation of one of 
these quantities with l m . We choose a = g(0), for 
which our theory predicts the leading asymptotic behav- 
ior g(0) = (In 2 l m )/(2w 2 l m ). In Fig. 7, we plot the nu- 
merically calculated values of g(0) versus l m and compare 
with the theoretical prediction. The overall shape of the 
Z rra -dependence is well captured by the theory, although 
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FIG. 6. Parametric plot of the nematic shadow density p N 
against p, for the same parent distribution and cutoff range 
as in Figs. 4 and 5. Very good agreement is observed between 
the numerical results (solid) and the theoretically predicted 
relation (dashed). 
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FIG. 7. Plot of numerically calculated values of g(0) 
against l m (solid) together with the leading asymptotic behav- 
ior g(0) = (In 2 l m )/(2w 2 l m ) predicted theoretically (dashed), 
for a log-normal distribution with a = 0.5. The agreement is 
satisfactory, though corrections to the asymptotic theory are 
clearly still important in the range of l m shown. 



sublcading corrections, which from Eq. (37) are of rel- 
ative order C(l/lnZ m ), are clearly not yet negligible in 
the range of l m considered. In summary, then, all numer- 
ical results are consistent with the theoretical predictions 
derived above. 



V. THE SCHULZ DISTRIBUTION 

Having observed the rather surprising effects caused 
by the long rods at the onset of isotropic-nematic coexis- 
tence in systems with fat-tailed length distributions, an 
obvious question is whether similar phenomena arc possi- 
ble even for more strongly decaying length distributions. 
We therefore now analyse, using the same numerical and 
theoretical methods as above, the case of a Schulz dis- 
tribution of lengths. For this distribution our previous 
studies of the Zwanzig 22 and V2 Onsager models 27 did 
not show any signs of the phase behavior being driven 
by the long rods. However, comparing our above results 
for the log-normal case with those obtained for the V2 
Onsager model 26 , it is clear that in the full Onsager the- 
ory the effect of the long rods is much more pronounced 
than in the approximate models. Long rod effects might 
therefore also appear, in the full Onsager theory, for the 
more strongly decaying Schulz distribution, but would 
then be expected to be weaker than for the log-normal 
case. 

The Schulz length distribution can be written as 



p(0)(0 = T(7Tiy Pcxph(z + 1)/] 



(42) 



where we have again imposed an average length of 1 (for 
infinite cutoff l m ). The polydispersity a, defined as be- 
fore as the relative standard deviation of the distribution, 
is related to z by 



(I 2 ) ~ {I? 
(I) 2 



1 



(43) 



From Eqs. (14) and (42) it is clear that if g(0) > z + 1 
the nematic density distribution is again exponentially 
diverging for large I. Assuming initially that this is not 
the case, however, one can solve numerically Eqs. (12) 
and (13) for the onset of isotropic-nematic phase coex- 
istence. The results for g(0) are shown in Fig. 8. For 
large z, i.e., small polydispersity a = (1 + z) -1 / 2 , one 
has g(0) < z + 1 (dashed line). In this regime the ne- 
matic density distribution decays exponentially for large 

1, and the results are essentially independent of the cutoff 
l m , which could in fact be taken to infinity. For smaller 

2, on the other hand, g(0) > z + 1, implying that the 
nematic density distribution is exponentially increasing 
and a finite cutoff l m is necessary. In this regime, the 




FIG. 8. The parameter g(0) of the nematic shadow phase, 
plotted against z for a Schulz distribution with cutoff 
l m = 100. See text for discussion. 
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situation resembles the case of the fat-tailed length dis- 
tributions discussed earlier, where for large enough cutoff 
and polydispersity the less than exponentially decaying 
length distribution was not able to balance the diver- 
gence of the factor exp[Z <?(())] in Eq. (14). With a Schulz 
distribution, the only difference is that now the compari- 
son is between two exponential terms exp[— (z + 1)1] and 
cxp[Z(7(0)]. Given this analogy, it is not surprising that 
the cloud and shadow curves (Fig. 9) show behavior qual- 
itatively similar to that found for the log-normal distri- 
bution: a kink in the cloud curve and a discontinuity in 
the shadow curve again indicate the presence of a three- 
phase I-N-N coexistence region in the phase diagram. 
Quantitatively, however, the kink in the cloud curve 




FIG. 9. (a) Cloud curves for Schulz distributions with cut- 
off l m = 50 (solid) and l m = 100 (dashed), plotted as cloud 
point density p versus polydispersity a on the y-axis. The 
system is again dominated by the long rods for large cr, but 
now the dependence on the cutoff is much less pronounced 
than for a log-normal length distribution (Fig. 1). The cloud 
curves also exhibit a kink, more clearly visible only on a mag- 
nified scale (inset). The dotted line represents the theoreti- 
cally predicted limiting form p = l/(ci a ) of the cloud point 
curve above the kink (see text) in the limit l m — > oo. (b) 
Corresponding shadow curves. The discontinuity correspond- 
ing to the kink in the cloud curves is much narrower than 
in the log-normal case (Fig. 1), and is visible clearly only in 
the inset. The dotted line again gives the theoretical predic- 
tion for the shadow curve above the discontinuity in the limit 
l m -> oo, p N = (1/a 2 + l/2a 4 )/(3 Cl ). 

(Fig. 9(a)) is now much less pronounced, and the cutoff- 



dependence of the cloud curve above the kink is also 
rather weaker. Similar comments apply to the shadow 
curve (Fig. 9(b)): the discontinuity is still present but 
very small, with the nematic phases on the two different 
branches having very similar densities (Fig. 9(b), inset). 
As for the log-normal distribution, for small polydisper- 
sities (i.e., below the kink or discontinuity, respectively) 
the cloud and shadow curves are essentially independent 
of the cutoff and connect smoothly with the monodis- 
perse limit at a = 0. 
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FIG. 10. Scaled volume fraction representation of the 
shadow curve for Schulz distributions with cutoff l m — 50 
(solid) and l m = 100 (dashed). Notice that, although the 
discontinuity is much more visible here than in the density 
representation (Fig. 9(b)), the maximum value of pi is not 
reached at the discontinuity. 

The discontinuity in the shadow curves is much more 
apparent in the scaled volume fraction representation 
(Fig. 10). Notice that now the discontinuity does not 
coincide with the point where pf reaches its maximum 
value, as was the case for the log-normal distribution 
(Fig. 2). This is also reflected in a plot of the average 
rod length in the nematic shadow against polydispersity 
(Fig. 11), which at the discontinuity in the shadow curve 
jumps to a large value but one that is still some way be- 
low l m . In this region the nematic phase thus contains 
many long rods, but is not yet entirely dominated by only 
the longest rods. This can also be understood by looking 
back at Fig. 8: for values of z just below the discontinuity 
(corresponding to the values of a just above the disconti- 
nuity in Figs. 10 and 11), ^(O) is very close to z + 1 and 
so the overall exponential factor exp{[g(0) — (z + l)]l]} 
is almost constant over the range I = . . . l m , giving a 
broad nematic length distribution (14) dominated by the 
non-exponential factors I* J d8 e 1 ^- 9 ^" 9 ^ . 

A. Theory for Schulz distributions with large cutoff 

From the numerical results obtained above, it looks 
like the case of the Schulz distribution is actually rather 
similar to the log-normal case. In both cases, the pres- 
ence of the long rods strongly affects the phase behavior 
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FIG. 11. Average rod length in the nematic shadow phase 
for Schulz distributions with cutoff l m = 50 (solid) and 
l m = 100 (dashed), plotted against polydispersity a on the 
y-axis. Above the discontinuity, the average length jumps to 
a large value but remains below l m , implying that the nematic 
phase is not yet entirely dominated by the longest rods. 



above a certain value of the polydispersity a; the thresh- 
old value of a tends to zero as l m increases for the log- 
normal case, but appears essentially independent of l m 
for the Schulz distribution. Above the threshold, the ne- 
matic phase is dominated by the long rods present in the 
system, although for the Schulz distribution the average 
length seems to remain rather below l m . Given these 
similarities, we now investigate whether the theory that 
we developed for the log-normal case can be extended to 
the case of the Schulz distribution. As we will see, the 
central assumption of dominance of the long rods in the 
nematic can still be made self-consistent. 

If long rods again dominate the density distribution 
in the nematic phase, we can repeat all the steps up to 
Eq. (25) in Sec. IV. The same equation for h{ff) follows, 
and if we assume that p c ff = p^l^ is large we get the 
scaling solution h(6) = h(p c gs'm9). Angular integrals 
can then again be simplified because only small values 
of 6 ~ 1/poff are relevant, leading to Eq. (30), to the 
simple result II N = 3p N for the osmotic pressure in the 
nematic phase, and to the expression (32) for g(0) = 
c\p — (8/n)(t) t /l m - At this point the different shape of 
the Schulz distribution enters: for large l m we now expect 
that 5(0) -> z + 1, rather than g(0) -> 0. This then 
implies that the cloud point density p has a finite limit 
for large cutoff, given by 
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(44) 



from Eq. (43). From the osmotic pressure equality 3p N = 
p+ (c 1 /2)p 2 the nematic shadow density then also has a 
finite limit, 
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(45) 



These theoretical predictions for the limiting form of the 
cloud and shadow curves in the cutoff-dominated regime 
are shown by the dotted lines in Fig. 9, and are certainly 
plausible given the numerical results for finite cutoffs. 

To obtain the leading terms in the approach of p and 
p N to their limiting values, and to establish the threshold 
value of polydispersity a above which the onset of I-N 
coexistence is affected by the presence of the long rods, 
let us define 5 = g(0) — (z + 1). As pointed out above, 
Eq. (30) still holds for the Schulz distribution case, but 
now yields to leading order 
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where we have used P^°'(l) cx l z e~^ z+1 ^ 1 and c collects 
all numerical constants as well as the factor p/(p N ) 2 
which approaches a constant for l m — > 00 from Eqs. (44) 
and (45). If d converges to zero slowly enough with 
l m — > 00 for the product Sl m to diverge, then the ex- 
ponential factor is dominant in the integral in Eq. (46) 
and we can replace l z by (see App. C) to get 
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Rearranging gives, with a new constant c' which contains 
the finite limit value of p N , 
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The inverse of the function of Sl m on the l.h.s. is asymp- 
totically just a logarithm, yielding 



S ~ (3 - z) 



In/ 



(47) 



For z < 3, corresponding to a > 1/2, this result is con- 
sistent with our assumptions: we have d — > and thus 
g(0) — > z + 1 as expected, and also Sl m — > 00 as assumed 
above. From the convergence of S to zero one can then 
also obtain the approach of p and p N to their limit val- 
ues: e.g. cip = g(0) + (8/TT){t) t /l m = z + l + S+0(l/l m ), 
giving p — {z + l)/ci ~ 5 to leading order. 

From the above results we can furthermore estimate 
the average rod length in the nematic shadow phase for 
z < 3 and finite l m . The exponential factor cxp(Sl) dom- 
inates the nematic density distribution and so rods in a 
range of 0(1/ 8) below l m should contribute to the aver- 
age length, giving the estimate 
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This shows that there are strong logarithmic corrections, 
consistent with the fact that in Fig. 11 the average ne- 
matic rod lengths (Z) N are still significantly below l m . 
By contrast, in the log-normal case, where the role of 
5 is played by g(0), the relative corrections to (Z) N are 
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~ l/(g(0)l m ) ~ l/ln 2 Z m ; even though still logarithmic, 
these correction terms are significantly smaller. 

Overall, the behavior of "wide" Schulz distributions 
with a > 1/2, i.e. z < 3, is therefore rather similar to that 
of log-normal distributions. We again have (Z) N — ► l m for 
large l m , although now the corrections are rather more 
important than in the previous case, and the rods in the 
ncmatic are strongly ordered (since p c ff = P N ^m ~ lf n 
diverges for l m — ► oo). The main difference is the fact 
that the cloud and shadow densities, p and p N , now tend 
to distinct nonzero limits for l m — ► oo rather than to zero: 
the smaller number of long rods in the Schulz distribution 
is not sufficient to induce phase separation at arbitrarily 
small densities. 

So far we have only covered the regime z < 3. The 
case z = 3 requires a more careful treatment; here the 
leading contribution to S calculated in Eq. (47) vanishes, 
and it turns out that S scales as l/l m , with 6l m ap- 
proaching a finite limit rather than diverging. Rods in a 
range ~ 1/6 ~ l m now contribute to the nematic density 
distribution, which therefore is no longer dominated by 
the longest rods alone even for l m — > oo; instead, one 
finds that the distribution approaches a scaling function 
of l/l m , with the ratio (/) N jl m approaching a nontrivial 
limit value < 1. 

In the case z > 3, finally, we cannot construct a self- 
consistent theory based on the assumption that the ne- 
matic phase is dominated by long rods. Looking at 
Eq. (46), one sees that for z > 3 negative values of 6 
would be required to make the r.h.s. of the equation (and 
thus p N ) finite for l m — > oo; with such negative values of 
5, the assumption of dominance of the long rods in the 
nematic phase is no longer self-consistent. One might try 
the milder assumption that the nematic is dominated by 
rods which are long but still short compared to l m , as- 
suming e.g. (Z) N ~ Z° with a < 1. We found, however, 
that this always leads to contradictions. Our conclusion 
is therefore that for z > 3, i.e., a < 1/2, the nematic 
phase is dominated by "short" rods with lengths not in- 
creasing with l m . The theoretically predicted threshold 
value of a below which the cutoff l m is irrelevant and the 
presence of long rods in the system does not significantly 
affect the phase behavior is therefore a = 1/2. 

Our numerical results (Fig. 9) show the kink in the 
cloud curve, and the corresponding discontinuity in the 
shadow curve at a m 0.48, i.e., close to but slightly below 
the theoretically predicted threshold value. Comparing 
with Fig. 8, this corresponds to the fact that, coming 
from large z, the discontinuous jump in <j(0) occurs before 
the extrapolation of the solution in the large z-regime 
intersects the critical line <?(0) = 2+1; consistent with our 
theory, and extrapolating by eye, this intersection occurs 
close to z — 3. It therefore appears that for finite l m the 
shadow curve jumps to a cutoff-dependent branch, i.e., 
a nematic containing many long rods, already slightly 
below the asymptotic threshold value a = 1/2. The value 
of cr where this jump occurs should then increase towards 
1/2 as l m increases. 



VI. CONCLUSION 

We have studied the effect of length polydispersity on 
the onset of I-N phase coexistence in the Onsager the- 
ory of hard rods. To assess the possible effects of long 
rods, two different length distributions were considered, 
one with a slowly decaying, "fat" tail (log-normal) and 
another with an exponentially decaying tail (Schulz). 

The log-normal distribution was chosen because it had 
earlier given interesting results in the context of ho- 
mopolymers with chain length polydispersity 25,24 and in 
our previous investigation of the V2 Onsager model 26 , 
obtained from a truncation of the angular dependence of 
the excluded volume of the Onsager theory We showed 
that a length-cutoff l m needs to be introduced for fat- 
tailed distributions such as the log-normal to avoid diver- 
gences in the equations for the onset of phase separation; 
the presence of such a cutoff is of course also physically 
reasonable. The most striking result from our numeri- 
cal solution for the properties of the isotropic cloud and 
ncmatic shadow phases is that the cloud curves show a 
kink and the shadow curves corresponding discontinu- 
ities: this establishes that for fat-tailed unimodal length 
distributions three-phase I-N-N coexistence occurs. The 
cloud and shadow curves show a strong dependence on 
the cutoff length, with both moving rapidly to lower den- 
sities as the cutoff increases. A plot of the average rod 
length in the nematic shadow phase suggested that it 
consists almost entirely of the longest rods in the sys- 
tem, i.e., those of length l m ; as a result, the nematic 
phase also exhibits very strong orientational ordering. 

A theoretical analysis of the limiting behavior for 
l m — ► 00 confirmed and extended these numerical re- 
sults. For large cutoffs, the ncmatic indeed comprises 
only the longest rods in the parent length distribution, 
and is very strongly ordered. Beyond this, the theory also 
predicts that the densities of the isotropic cloud and ne- 
matic shadow phases in fact vanish (with constant ratio 
p/p N = 3) in the limit of infinite cutoff. This rather 
surprising result means that even though the average 
rod length in the parent distribution is finite, the fat 
tail of the distribution ensures that enough arbitrarily 
long rods are present to induce phase separation at any 
nonzero density. Even though the ncmatic shadow den- 
sity converges to zero for increasing l m , it does so slowly 
enough (p N ~ (In 2 l m )/l m ) for the rescaled rod volume 
fraction of the nematic to diverge logarithmically with l m 
(pf ~ In 2 l m ). For any given aspect ratio of rods, the On- 
sager second virial approximation thus eventually breaks 
down as l m increases, at the point where the true vol- 
ume fraction (D/Lq)p^ of the nematic shadow becomes 
non-negligible compared to unity. 

We then studied the case of a Schulz distribution of 
rod lengths, for which our previous studies of the sim- 
plified Zwanzig 22 and V2 Onsager models 27 showed no I- 
N-N coexistence and no unusual behavior in the limit of 
infinite rod length cutoff. The full Onsager theory stud- 
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ied here revealed, however, that such effects do indeed 
occur: above a threshold value of the polydispcrsity a, 
the numerical results show that the nematic density dis- 
tribution becomes exponentially divergent for large rod 
lengths; a finite cutoff l m again needs to be imposed to 
get meaningful results. Above the threshold, the cloud 
and shadow curves then depend on l m , although much 
more weakly than in the log-normal case. The dominance 
of the long rods in the nematic shadow phase above the 
threshold is also weaker than for the log-normal, with 
average rod lengths that are large but significantly be- 
low l m . At the threshold itself, a kink in the cloud curve 
and a discontinuity in the shadow curve occur, indicat- 
ing that even for the "well-behaved" Schulz distribution 
the Onsager theory predicts a three-phase I-N-N phase 
coexistence region in the phase diagram. 

We were again able to clarify these results by theoret- 
ical analysis of the limit l m — > oo. We found that the 
limiting threshold value of the polydispersity is a = 1/2, 
corresponding to an exponent z = 3 in the Schulz dis- 
tribution, in good agreement with the numerically cal- 
culated thresholds at small cutoffs. Above the threshold 
value theory predicts that the average nematic rod length 
approaches l m as in the log-normal case, but now with 
larger logarithmic corrections which explain the smaller 
average lengths observed numerically. In contrast to 
the log-normal distribution the cloud and shadow curves 
above the threshold approach finite limiting values for 
l m — > oo. The physical interpretation of this result is 
that the smaller number of long rods in the Schulz dis- 
tribution is not sufficient to induce phase separation at 
arbitrarily small densities. 

It is appropriate at this stage to compare the above re- 
sults with our earlier analysis of the T>2 Onsager model 26 . 
For the Schulz distribution we have mentioned already 
that the Vi Onsager model, in contrast to the full On- 
sager theory, predicts no unusual effects (I-N-N coexis- 
tence and cutoff dependences) due to the presence of long 
rods. For the log- normal distribution, the Vi Onsager 
model does exhibit a three-phase I-N-N coexistence, with 
cloud curves showing a kink and shadow curves a corre- 
sponding discontinuity. As in the full Onsager theory, 
above the kink/discontinuity the cloud curves are also 
strongly dependent on the cutoff value. However, the 
limiting behavior of cloud and shadow curves is rather 
different: both the densities and the rescaled rod vol- 
ume fractions of the isotropic cloud and nematic shadow 
phases converge to finite, and in fact identical, limiting 
values for large cutoff. The nematic was also not dom- 
inated by the longest rods. In fact, the isotropic and 
nematic phases differed only through a somewhat larger 
fraction of long rods contained in the nematic, and the 
shorter rods show only negligible order in the nematic, 
with the overall orientational order parameter vanishing 
in the limit. The enhancement of long rods in the ne- 
matic only becomes detectable in higher order moments 
of the nematic density distribution such as J dl /™p N (7), 
which diverge for n > 3/2. 



The above differences between the predictions of the Vi 
Onsager model and the full Onsager theory can be un- 
derstood as follows. In the Onsager theory, the excluded 
volume of two rods vanishes as the angle between the 
rods decreases to zero; this favors strongly ordered ne- 
matics such as the nematic shadow phase dominated by 
long rods that we found at the onset of phase coexistence. 
In the Vi Onsager model, on the other hand, and indeed 
in any similar truncation of the expansion of the kernel 
K(0, 6') in Lcgcndre polynomials 27 , the excluded volume 
remains nonzero even for two rods fully aligned with the 
nematic axis. This disfavors nematic phases containing a 
substantial number of long and strongly ordered rods. It 
thus makes sense that the nematic shadow phase even for 
the "fat" log-normal distribution is predicted to contain 
only a small (though enhanced compared to the isotropic 
phase) fraction of long rods. 

Looking back over our results for the effects of length- 
polydispersity in the full Onsager theory, it is clear that 
all the effects of long rods that we observe arise from the 
exponential factor exp[7<?(0)] (see Eq. (14)) which dom- 
inates the enhancement of the nematic shadow phase's 
density distribution over that of the isotropic parent 
phase. Any parent length distribution with a less than 
exponentially decaying "fat" tail will therefore exhibit 
divergences in the nematic distribution, leading to phase 
behavior similar to that for the log-normal case. In fact, 
our theory in Sec. IV applies to all such fat-tailed dis- 
tributions. The Schulz distribution with its exponential 
tail is the borderline case, where one cannot predict a 
priori whether the presence of long rods will have sig- 
nificant effects. We found that it does, above a thresh- 
old value of the polydispersity. To our knowledge this 
is the first time that such an effect has been observed 
in polydisperse phase equilibria. In Flory-Huggins the- 
ory for homopolymers with chain length polydispersity, 
for example, where the enhancement factor is also a lin- 
ear exponential (in chain length), no long rod-effects are 
found for Schulz distributions 24,25 . Finally, for parent 
rod length distributions decaying more than exponen- 
tially, e.g. as ~ exp(— l a ) with a > 1, no cutoff de- 
pendences are expected since the nematic density dis- 
tribution will always be well-behaved for large lengths. 
Of course, this does not mean that I-N-N phase coex- 
istence is excluded for such distributions. Consider for 
example a log-normal length distribution modulated by 
a Gaussian factor exp[— ^ 2 /(2^)] with large l m . As we 
just saw, there is then no need for an explicit cutoff. On 
the other hand, the Gaussian factor will act as an ef- 
fective "soft" cutoff (hence the notation l m ). For large 
enough l m one thus expects phase behavior qualitatively 
similar to that discussed above for a "hard" cutoff, in- 
cluding I-N-N phase coexistence signalled by a kink in 
the cloud curve. 

Above, we have focussed exclusively on the onset of 
phase coexistence; both numerically and theoretically the 
analysis of the phase behavior inside the coexistence re- 
gion would be far more challenging. One question which 
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one would like to answer, for example, concerns the over- 
all phase diagram topology. On the one hand, the three- 
phase I-N-N region could be confined to a narrow density 
range inside the I-N coexistence region, as is the case for 
the Vi Onsager model with a log-normal length distri- 
bution 26 . The alternative would be for the I-N-N region 
to extend all the way across the I-N coexistence region, 
connecting to the ncmatic cloud curve and being bor- 
dered by a region of N-N coexistence; this is the behavior 
predicted by Onsager theory for bidispcrse systems 16,23 . 
Apart from a direct numerical attack on the phase coex- 
istence region, which for now seems out of reach, clues 
to an answer could be provided by an approach based 
on the Flory lattice model of hard rods 31 . For the sce- 
narios studied in the past this has yielded results qualita- 
tively similar to the Onsager theory 32 ' 33,15 , in spite of the 
rather crude treatment of the orientational entropy. Pre- 
liminary work shows that, the limit of thin rods, Flory's 
excess free energy corresponds to an excluded volume 
term which correctly tends to zero for small rod angles. 
Together with the full expression for the ideal part of the 
free energy this can be shown to give scaling behavior 
in the limit of strong ordering (i.e., high density) very 
similar to that of the full Onsager theory. This version 
of the Flory lattice model may therefore produce predic- 
tions that are more in qualitative accord with Onsager 
theory than e.g. the Vi Onsager model. It shares with 
the latter the desirable feature of being "truncatable" , 
having an excess free energy that depends only on two 
moments of the density distribution. This will allow ef- 
ficient calculation of phase equilibria using the moment 
free energy method 22 ' 27,34 ™ 37 ; work in this direction is in 
progress. 



APPENDIX A: HIGH-DENSITY SCALING FOR 
THE MONODISPERSE ONSAGER THEORY 

We summarize here the arguments leading to the 
scaling solution for the angular distribution in ncmatic 
phases at high density, for the case of monodisperse 
rods 38 . The relevant free energy is obtained from the 
polydisperse version (1) by dropping all ^-integrations 
and setting 1 = 1, giving 



f = p(lnp-l) + p J d6 P{6)\tlP{6) 
+ -p 2 (dOdl)' P(6)P(6')K(6,6') 



(Al) 



This expression needs to be minimized with respect to 
P( 9), subject to the normalization condition J d9 P(9) = 
1 in our usual notation. One obtains 



P{6) = 



- ip(6) = -p / d6> P(6')K(6,6') 
Jd6> eW) 1 



which is the obvious monodisperse version of Eq. (2). 
Defining the function h( 9) — ?/>(0) — ip(9), which obeys 
h(0) =0, this can be written as 



P(0) = 



,-h(9) 



(A3) 
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h{9) =p I dff P{9') [K(0, 6') - K(0, 9')} (A4) 



Now consider the regime of high densities, p 3> 1. From 
Eq. (A4), it is clear that for large density, h{9) becomes 
large. P(9) then becomes strongly peaked around 8 = 
so that the only non- vanishing contribution to the angu- 
lar integral in Eq. (A4) comes from the range 9' « 1. 
(Here and in the following we use the symmetry of P(6) 
and h(9) under 9 — > it — 8 to restrict all integrations to 
the range 6 = 0... it/2.) For 9 ~ 0(1) > 9', we can then 
approximate K(9, 9') ~ K(6, 0) = (8/tt) sin6> so that 



h(9) ~ — psinf 

TT 



(A5) 



to leading order in p. This expression will not be valid 
for small 9, but suggests that in this regime a scaling 
solution in terms of the scaling variable t = p sin 9 could 
exist, i.e. h{9) = h(t). For consistency with Eq. (A5) for 
9 ~ 0(1) (and large p), the scaling function should then 
have the leading asymptotic behavior h(t) = (8/7r)i for 
large t. Now, written in terms of h(i), Eq. (A4) reads 

h(t)= ^-W'F = (A6) 

r dt' v e -h(t>) 

where the integrals are over the range . . . p and 9 = 
arcsin(i/p), 9' = arcsin(t'/ j o). The key property that 
allows one to get a density-independent equation for h(t) 
is the scaling behavior of the kernel. For finite t and t' 
and large p, one has 9 m t/p, 9' m t'/p and for such 
small (and comparable) angles the kernel scales linearly 
with the angles. The product pK(9, 9') thus approaches 
a finite limit for p — > oo, 

pK(t/p, t/p') - K(t, t') = V^+^F , (A7) 
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In the same limit we can replace the factors [1 — 
(t'/p) 2 ] -1 / 2 in Eq. (A6) by 1, and obtain the scaling equa- 
tion 



hit) = {K{t,t') - K(0,t%, 



(A8) 



where the average over t' is over the normalized proba- 
bility distribution 



(A2) 



P(t) = -te~ h ^\ 7 - / dtte- h ® 
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of the scaling variable t, now running over the range 
. . . oo since we have taken p — * oo . A plot of the nu- 
merical solution of Eq. (A8) is shown in Fig. 12 and has 
the expected leading linear behavior at large t. 



with large cutoff. The excess free energy is the last term 
in Eq. (All), thus 



-p{K(t,t%, t ,=2p 



(A13) 




FIG. 12. Numerical solution of Eq. (A8) for the scaling 
function h(t) (solid) together with the asymptotic linear be- 
havior at large t (dashed). 



In the high-density limit, by the same arguments as 
above, the normalization factor for the orientational dis- 
tribution (A3) becomes J d6 exp[— h(psm6)] = 7/p 2 , 
with 7 the normalization factor defined in Eq. (A9). Thus 
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Inserting this into the expression (Al) for the free energy 
and transforming everywhere from 9 to t, one has 

f = p{\np-l)+2p\np 

- p [(h(t)) t - ln 7 ] + \p{K{t, t')) t , t , (All) 

From this it follows that the osmotic pressure for large 
densities is simply 



n 



df , _ 
P dp- f = 3p 



(A12) 



It then also follows that the excess free energy / (the 
last term in Eq. (All)) is just 2p. This can be seen by 
comparing the result (A12) with that obtained via a dif- 
ferent route: since P(9) is determined to minimize the 
free energy, one can evaluate II = p df /dp — f by differ- 
entiating Eq. (Al) while holding P(9) constant. Because 
the excess free energy is quadratic in p for constant P{&), 
this gives n = p + f and therefore / = 2p by comparison 
with Eq. (A12), as claimed. 

The result / = 2p is also useful for deriving an identity 
which we use in the main text to show that II = 3p holds 
in the nematic shadow phase at the onset of phase sep- 
aration in a system with a fat-tailed length distribution 



From Eq. (A8) and K(0,t') = (8/tt)*' we also have 
(K(t,t')) t/ = (K(t,f) - K(0,t%, + ~{t') v 

TT 

= h(t) + -(t) t (A14) 

Inserting into Eq. (A13) we then obtain the desired iden- 
tity 



(M*)>* + "(*)*= 4 

TT 



(A15) 



Notice that arguments very similar to those above ap- 
ply also to polydisperse nematics: one again has a scal- 
ing solution h{9) = h(t) for high density, in terms of the 
same scaling variable, and this again leads to II = 3p and 
/ = 2p. Van Roij and Mulder 18 showed this explicitly for 
the bidisperse case. 



APPENDIX B: THE APPROXIMATION OF 
DOMINANCE OF LONG RODS 

Our theory in Sec. IV for the onset of phase separa- 
tion in systems with fat-tailed length distributions was 
based on the assumption that the nematic shadow phase 
is dominated by the longest rods in the system. This al- 
lowed us to replace terms that depended weakly on rod 
length I by their values at the cutoff l m . We now verify 
that this assumption is justified in the four cases where 
we have used it, namely in Eqs. (21), (22), (23) and (25). 
Let us start from the simplest of these, Eq. (21). Define 
the angular integral 



A(l) = / Me-W^W) 



The density distribution (20) in the nematic shadow 
phase is then p N (l) = pP(°) (l)e lg W A(l) . The normal- 
ized nematic length distribution P N (l) = p N (Z)/p N can 
be written as P N (Z) = Q(l)A(l)/A(l m ) if we define 



Q(l) 
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) 



Looking back at Eq. (21), the assumption of long rod 
dominance amounted to replacing the weakly varying 
factor A{l)/A{l m ) by A(l m ) / A{l m ) = 1, effectively sub- 
stituting Q(l) for P N (l). To check that this is justified, 
we need to consider the unapproximatcd Eq. (20), which 
after integration over I and division by p N reads 
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dl Q(l) 
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We effectively approximated this by J dl Q(l) = 1, so we 
need to show that the contribution from the short rods 
to Eq. (Bl) is negligible compared to unity. We will need 
the /-dependence of A{1) to do this. Restricting the inte- 
gration range to . . . 7r/2 by symmetry, using the scaling 
form of h(9) = h(t) and transforming to the scaling vari- 
able t = poft sin 6 gives 



A{1) 



PeS 



(It 



t 



-(l/l m )h(t) 



V 1 - WPeftY 



(B2) 



In the limit / — > one has A(l) = 1, and A(l) will remain 
of this order while the exponential exp[— (l/l m )h(t)} is 
close to unity even for t = p e g. Since p e ff is large and 
h(t) linear in t for large arguments, this gives the crite- 
rion (l/l m )PeS ~ 1, or I ~ Im/PcS- Up to this value of 
/, we can approximate A(l) ss 1. For larger /, the fac- 
tor exp[— (l/l m )h(p e ff)} is small enough for the integral 
to be dominated by values t <C p e ff, so that we can set 
[1 — (t/ pcfj) 2 ] 1 / 2 ss 1 and extend the upper limit of the 
^-integral to infinity. The bulk of the integral still comes 
from values of t 1, however, where h(t) is linear, and 
carrying out the integral (B2) with this approximation 
gives the scaling A(l) ~ p~g (/ m //) 2 . As / increases, the 
range of t- values contributing to the integral reduces, and 
eventually the quadratic behavior of h(t) near t = leads 
to corrections to this scaling. Even at I = l m , however, 
these effects are relatively small since h(t) is approxi- 
mately linear even down to t « 1 (see Fig. 12); we there- 
fore neglect them to a first approximation. In summary, 
we thus have the scaling A(l) ~ 1 for / <C l m /p c «, and 
A(l) — p~£{l m /l) 2 for / > l m /PcS- A sample plot of A(l) 
evaluated numerically together with our approximation 
is given in Fig. 13. 




le-05 0.0001 0.001 0.01 0.1 



10 100 



FIG. 13. Plot of A(l) against / for l m = 500 and 
p c ff = 6500 (~ / m ln 2 / m ). Notice the Jp^gHm/l) 2 scaling at 
intermediate values of I and the quasi-constant behavior at 
small I. At I ~ l m corrections to the (l m /l) 2 behavior are 
visible. 

It follows from the above results that the factor 
A(l)/A(l m ) in Eq. (Bl) is no larger than ~ (l m /l) 2 , 
even for very small /. The contribution to the /-integral 



from rod lengths / of order unity is therefore bounded 
by ~ Jdl Ifj- 2 Q{1). For I = 0(1) we can set the factor 
cxp[/g(0)] in Q(l) to one since g(0) is small (for large l m ); 
the factor p/p N is also asymptotically just an unimport- 
nat constant. The short rods contribution to Eq. (Bl) is 
thus of order ~ A(l m )lf n J dl P^{l)l~ 2 or, if we extend 
the integral up to / = oo, A^m)^ (l~ 2 ) with the average 
taken over the parent distribution. Now A(l m ) ~ 
(more precisely, A(l m ) = r yp~g] see before Eq. (AlO)) 
so this contribution scales as (/ m /Poff) 2 ~ {ImP^) 2 - For 
the log-normal length distribution, p N ~ (In 2 l m )/l m , and 
so the short rods contribution ~ l/ln 4 /„,, to Eq. (Bl) 
does indeed become negligible for large l m , when com- 
pared to the long rods contribution of 1. 

Next consider the approximated osmotic pressure 
equation (25). We need to show that the contribution 
of the short rods to the /-integral in Eq. (24) is negligible 
compared to the long rods part, which we evaluated to 
be 2,o N . Dividing by p N to have a quantity to which the 
long rods contribute a value of order unity, and discard- 
ing factors which are constant for large l m , we have to 
consider the integral 

dl W^^C)-^] 



A(l) 



The first term in the curly brackets is easy: the /-integral 
is proportional to / dl Q(l)(l/l m )[A(l)/A(l m )]. Compar- 
ing with the integral J dl Q(l)[A(l)/A(l m )] treated above, 
the short rods contribution here is suppressed by an ad- 
ditional factor of l/l m and so definitely negligible. The 
second term, on the other hand, is of the form 



dl Q(l)l 



A{1) A' (I) 
A(l m ) A(l) 



dl Q(l)l 



A(l) d 
A{l m ) dl 



In A(l) 



As we saw above, A(l) varies at most as a power-law with 
/, so that (d/dl)lnA(l) ~ 1// and we are again led back 
to the integral / dl Q(l)[A(l)/A(l m )] for which we showed 
the dominance of the long rods above. (This argument 
applies even in the small-/ range / < l m / p B s, where A(l) is 
approximately constant and (d/dl) In A(l) therefore even 
smaller.) 

In Eq. (18), which we approximated by Eq. (22), we 
have a similar /-integral that turned out to scale as l// m 
(see Eq. (32)). Multiplying then by /,„ to again have a 
long rods contribution of order unity, and using pll m = 
(p/p N )(l/lm)PcB ~ {l/lm) Pea, wc nc °d to consider the 
integral 



dl Q(l) 



I All) 

'lit 

A(l 



d6 — p c ff sin 9 

7T 



-(l/l m )h(p tt sin 9) 



A{1) 



(B3) 



As before, the angular integral will be dominated (with 
the exception of very small lengths / < / m /p e ff ; see below) 
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by the range where t = p e g sin 9 is large. In this range 
h{t) as {8/ir)t is linear to a good approximation, so that 
the angular integral can be written as 



d9 h(p e g sin 9)- 



-{l/l m M PeS sine) 

Ml) 



A(l) 



The overall integral (B3) thus be- 
comes / dl Q(l)l[A(l)/A(l m )](d/dl)\nA(l) for which long 
rods dominance has been shown already. The contribu- 
tion from very small l-values I < l m j p e ff needs to be 
treated separately: here we use the fact that the angu- 
lar integral in Eq. (B3) can be viewed as an average of 
(8/ir)p e s sin# over a normalized distribution, giving a re- 
sult of at most ~ p c ff. Using that also A(l) ~ 1 in this 
regime, one needs to integrate p e aQ(l)(l/l m )[^/ A(l m )\ 
over the range I = 0...l m / p e g; the integrand is bounded 
there by Q(l)/A(l m ) ~ P(°)(Z) and gives a vanishing 
integral since the upper limit of the integration range, 
l m / Peg = l/(p N l m ) ~ l/\n 2 l m vanishes in the limit 
l m — > oo. 

Finally, we need to analyse Eq. (19), which we approx- 
imated by Eq. (23). Let us rewrite the r.h.s. of Eq. (19) 



Peg 

A(l m ) 



de>[K(6,e')-K(0,6')} f dlQ(l)^-t 



-(i/i m )h(e') 



One can view this as an average of the term in square 
brackets over an (unnormalized) distribution over 9' . We 
thus need to show that the short rods contribution to 



dl Q(l) 



-(l/l m )h(8') 



(B4) 



is negligible compared to the long rods contribution, at 
least for the values of 9' that are in the bulk of this 
distribution (rather than the tail). When evaluating 
the long rods contribution it is not a priori clear that 
one can treat the exponential factor exp[— (l/l m )h(9')] 
as weakly varying with I. We will see below that this 
can nevertheless be justified, so that the long rods con- 
tribution is simply e~ ht - e - 1 as used in Eq. (23). The 
short rods contribution to Eq. (B4) can be estimated 
by again approximating exp[/p(0)] as 1 and bounding 
cxp[— (l/l m )h{9')] < 1. This yields a contribution of 
at most / dl p( 0) (l)(l/l m )A(l m ) < (l)A(l m )/l m with the 
average again taken over P(°'(l) and thus giving unity. 
This is comparable with the long rods term exp[— h{9')\ 
only for values of 9' such that e~ h ^ e ^ ~ A(l m )/l m ~ 
~ C('m 3 m ~ 4 lm)- This means that the cor- 
rections due to the short rod integral become impor- 
tant only where the angular distribution has already de- 
cayed to ~ 1/lfn (up to logarithmic terms) of its value at 
9' = 0, i.e., far in the tails of the ^'-distribution. We can 
now also justify neglecting the /-dependence of the term 
exp[— (I / l m )h{9')\: at the onset of corrections from the 
short rods one has h(9') ~ lnl m . Up to this point, i.e., 
in the bulk of the distribution, the ratio h(0')/l m is at 



most ~ (\nl m )/l m . In the overall exponential factor in 
Eq. (B4), exp{/[<?(0) — h(9')/l rn ]}, this term is therefore 
still negligible compared to g(0) ~ (In l m )/l m and the 
resulting /-dependence can be ignored. 



APPENDIX C: THE DOMINANCE OF LONG 
RODS FOR THE SCHULZ DISTRIBUTION 

In this appendix we outline briefly how, for the Schulz 
distribution (42), one can establish that the assumption 
of the long rods being dominant in the nematic shadow 
phase is again justified (for z < 3). By way of example 
we discuss Eq. (21), which is obtained from Eq. (20) if 
long rods dominate. 

As we saw in App. B, the scaling of the angular integral 
A(l) =Jd9e-^ l ^ h W is p-£(l m /l) 2 for I > l m / PcS , and 
A(l) as 1 for smaller I. Inserting this into the nematic 
density distribution (20), p N {l) = pP^°\l)e l9 ^A(l), us- 
ing p e ff = p 1^ and exploiting that p and p N approach 
constant limits for l m — ► oo, one has p N (0 ~ l^l z ~ 2 e si 
for I > l m / Pes ~ an d p N (0 ~ l z e si for smaller I. 

We want to show again that the integral of p N (l) over 
the short rods (with lengths I of order unity) is negligi- 
ble compared to the long rods contribution (which, since 



P 



0(1), is of order unity). In the short rods regime we 



can approximate e si ~ 1. The integration of p N (l) over 
the range . . . l/l m then gives just ~ C(Z~ 2_1 ) which is 
negligible compared to unity for z > 0. The contribution 
from the range I > l/l m can be bounded by extending 
the integration not just over the short rods but in fact 
up to l m , giving ~ l^(l z ~ x - l x ~ z ) = {If' 3 ~ l^n' 1 ) 
which is again negligible compared to unity as long as 
z < 3. (For z = 1 the integral has a logarithmic cor- 
rection, giving ~ Z~ 2 lnZ m , but is still negligible.) As in 
the log-normal case, the integral over the nematic density 
distribution is therefore dominated by the longest rods, 
justifying Eq. (21). 
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